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Within a relativistic hydrodynamic framework, we use four different equations of state 
of nuclear matter to compare to experimental spectra from CERN/SPS experiments NA44 
and NA49. Freeze-out hypersurfaces and Bose-Einstein correlation functions for identical 
pion pairs are discussed. We find that two-pion Bose-Einstein interferometry measures the 
relationship between the temperature and the energy density in the equation of state during 
the late hadronic stage of the fireball expansion. Little sensitivity of the light-hadron data 
to a quark-gluon plasma phase-transition is seen. 
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The equation of state (EOS) of nuclear matter at very 
high energy densities remains unknown. A possible new 
state of nuclear matter, the quark-gluon plasma (QGP), 
may be formed within a very hot and dense zone of nu- 
clear matter - the fireball - during relativistic heavy-ion 
collision experiments. Many observables have been pro- 
posed as a signature for the QGP, among them mea- 
surements derived from intensity interferometry of iden- 
tical hadrons, also known as Bose-Einstein correlations 
(BEC). BEC functions are sensitive to the space-time 
dynamics of the fireball, and therefore should give clues 
about the EOS, which governs the evolution of those fire- 
balls. 

Among the many studies on BEC (for a recent overview, 
c/., the book by Weiner ]5J]) only few (c/., e.g., Rischke 
et al. [|| ) address the role of the EOS of nuclear matter 
in the intensity interferometry of identical hadrons. It 
is the purpose of this paper to further investigate the 
interplay between the EOS and BEC. In the following, 
we shall consider a framework of analysis which is based 
on relativistic hydrodynamics, since this approach allows 
for an explicit use of an EOS. 

We shall use the simulation code HYLANDER-C § with 
various equations of state, and Cooper-Frye M freeze- 
out. The results of these hydrodynamical (one- fluid- 
type) calculations will be compared to experimen- 
tal single particle momentum distributions from 158 A 
GeV Pb+Pb collisions, measured by the NA44 ||,|7| and 
NA49 @ Collaborations. After a discussion of the space- 
time features of the particular fireballs we shall compare 
BEC of identical pion pairs |9j to experimental correla- 



tion functions, as measured by NA44 . 
The first equation of state 111]], EOS-I, which we use in 
this study exhibits a phase transition to a quark-gluon 
plasma at a critical temperature T c — 200 MeV with a 
critical energy density e c — 3.2 GeV/ fm 3 . The second 
equation of state |lj], EOS-II, is a lattice QCD-based 
EOS (as is EOS-I) which has recently become very pop- 
ular in the field of relativistic heavy-ion physics. This 
equation of state includes a phase transition to a quark- 
qluon plasma at T c = 160 MeV with a critical energy 
density e c w 1.5 GeV/fm 3 . The third equation of state, 
EOS-III, has been extracted from the microscopic trans- 
port model RQMD jl3| under the assumption of com- 
plete thermalization, and does not include a transition 
to a QGP. We obtain a fourth equation of state, EOS-Ib, 
by changing the relationship between e c and T c in EOS-I 
to T c (e c = 1.35 GeV/fm 3 ) = 200 MeV. 
In Fig. 1 the four equations of state are shown in two dif- 
ferent representations. The particular parametrizations 
for the EOS which have been used here are described in 
Ref. H in more detail. 

In the following, we shall discuss five scenarios: we com- 
pare four calculations using EOS-I, EOS-II, EOS-III, and 
EOS-Ib, all at fixed freeze-out temperature Tf = 139 
MeV, and one calculation using EOS-I for fixed freeze- 
out temperature Tf = 116 MeV. 

It is possible to find initial distributions (c/. Table II 
in Ref. Q and Table I (here)) for the four equations 
of state, such that one can reproduce the single inclu- 
sive momentum spectra of 158 AGeV Pb+Pb collisions. 
For a freeze-out temperature Tf — 139 MeV, the initial 
conditions and a large number of final single inclusive 
momentum distributions for various hadron species have 
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been shown in Refs. 
sured by the NA44 



in comparison to the data mea- 
and NA49 [§ Collaborations. 
Those results refer to calculations using EOS-I, EOS-II, 
and EOS-III. Results showing the fits of single inclusive 
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particle momentum spectra using EOS-Ib with Tf = 139 
MeV and EOS-I with T f = 116 MeV and T f = 139 MeV 
are shown in Fig. 2. 

The magnitudes of the slopes in the transverse mass spec- 
tra have their origin in the freezeout temperature and 
the transverse velocity fields at freeze-out. Because the 
effective EOS softness || is larger in the calculations us- 
ing EOS-Ib with Tf = 139 MeV and EOS-I with T f = 
116 MeV compared to the calculation using EOS-I with 
Tf = 139 MeV, the values for the maximum transverse 
velocity at freeze-out, v™ ax , are correspondingly larger 
(cf. Table I) in the calculations using EOS-Ib with Tf = 
139 MeV and EOS-I with T f = 116 MeV compared to 
the calculation using EOS-I with Tf = 139 MeV. But 
because in the calculation using EOS-I with Tf = 116 
MeV the temperature is decreased while increasing the 
value for v™ ax , hardly any change is seen in the slopes of 
the transverse mass spectra compared to the calculation 
using EOS-I with Tf — 139 MeV. A more pronounced 
change is seen only in the calculation using the harder 
EOS-Ib while keeping Tf unchanged. 
It should be stressed that all calculations discussed so 
far (except the calculation using EOS-II) result in single 
particle momentum distributions that describe the data 
equally well. Although EOS-II was found in the calcula- 
tions of hadronic transverse mass spectra to be too soft 
(cf. Refs. HD), we shall use it here also for the calcula- 
tion of Bose-Einstein correlation functions. 
Before we discuss BEC, we briefly discuss the thermal 
evolution of the various fireballs. 

Fig. 3 shows the isothermes for the relativistic Pb+Pb 
fluids governed by EOS-I and EOS-Ib (cf. also Ref. §) 
until freeze-out has been reached. In Ref. Q it was shown 
that the calculation using EOS-I with T f = 139 MeV 
leads to a fireball of a much shorter liftime than the cal- 
culations using EOS-II and EOS-III with T f = 139 MeV. 
This behavior is caused by much smaller freeze-out en- 
ergy densities, €f, in the calculations using EOS-II and 
EOS-III compared to the calculation using EOS-I. We 
have e } = 0.292 [0.126 (0.130)] GeV/fm 3 when using 
EOS-I [EOS-II (EOS-III)]. A fluid that undergoes adi- 
abatic expansion needs more time to reach the smaller 
freeze-out energy densities. 

Since EOS-II and EOS-III yield similar lifetimes of the 
fireball, we attempt in the following to increase the life- 
time of the system which is governed by EOS-I. This can 
be achieved by (a) using a smaller freeze-out temperature 
Tf = 116 MeV, or (b) by hardening the EOS, i.e., using 
EOS-Ib instead of EOS-I (without changing T f = 139 
MeV). For the latter two cases, we obtain e/ = 0.127 
GeV/fm 3 (cf. Fig. 1 (c) and Ref. §). 
The lifetime of the fireball is reflected in Bose-Einstein 
correlations of identical pion pairs [O. Using the 



verse expansion |T(| as well as the contributions of reso- 
nance decay [[lTj in BEC. 

Fig. 4 shows data points taken by the NA44 Collab- 



oration 1 10 1, along with projections of calculated BEC 



functions for tt + 7t + pairs, using the same acceptance as 
the experiment. Consistent with expectation, the calcu- 
lations using EOS-II, EOS-III, and EOS-Ib give similar 
lifetimes and therefore sufficiently large longitudinally ex- 
panded fireballs (see Figs. 2,3). These calculations result 
in a good reproduction of the pionic NA44 BEC data. In 
addition, EOS-III and EOS-Ib produce an excellent de- 
scription |||5]|^] of hadronic single inclusive momentum 
spectra. 

It should be stressed here that a freeze-out temperature 
Tf = 139 MeV was adequate to achieve this agreement. 
Also consistent with expectation is that the calculation 
using EOS-I with Tj — 139 MeV gave a longitudinally 
expanded fireball which was too small. On the con- 
trary, it is initially surprising that a reduction of the 
freeze-out temperature to Tf — 116 MeV does not lead 
to a large enough longitudinal extension of the fireball 
in the calculation which uses EOS-I. The reason for this 
result is the following: a freeze-out temperature reduc- 
tion leads to a larger lifetime of the direct fireball, but 
because of the lower temperature, the relative fraction 
of heavy resonance decay contributions is reduced (by 
about 30%), so that the 'resonance halo' is reduced in 
size. Hence, the apparent fireball, which is a superposi- 
tion of the direct (or thermal) fireball and the resonance 
halo remains more or less unchanged in size. 
We note that the transverse projections of the correlation 
functions C2(q ut) and (^((jWe) are also described rea- 
sonably well, especially by those calculations which yield 
the larger longitudinal freeze-out extension. The num- 
bers in Fig. 4 can be obtained from a 1-dimensional fit 
of the correlation functions with jl8| 
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1 + \(K)-ex V [-<l?Rl(K) - q 2 R 2 (K) 
- q 2 s R 2 s (K) + 2 Pol (K)q oqi R (K)Ri(K)] 
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Bertsch-Pratt parametrization |15|, a sensitive quantity 
is the longitudinal projection of the two-pion correlation 
function, C2(qiong), since it reflects the effects of trans- 



but they should not be taken too seriously. In eq. (|l|) 
the qi (i = I, s,o) refer to the components of the momen- 
tum difference q = k% — k 2 , and K = ^(ki + k 2 ) is the 

average pion pair momentum. Furthermore, X(K) is the 
momentum dependent incoherence factor which accounts 
for reductions of the BEC due to long-lived resonances 
]l7j and averaging due to phase-space, respectively. 
In conclusion, by inspecting Fig. 1(b) we can see that 
only those equations of state which go through point 
C in Fig. 1 reproduce the experimental data on Bose- 
Einstein correlations fairly well. Changing the freeze- 
out temperature shows hardly any effect. Therefore, the 
measurements of Bose-Einstein correlations tell us which 
relationship between temperature and energy density is 
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neccessary for a valid choice of an equation of state in 
the calculations. Unfortunately, from the above consid- 
erations it must be noted that a two-particle BEC used 
by itself cannot be used as a tool to determine a possible 
phase-transition to a QGP, because the BEC show little 
sensitivity to the structure of the EOS. 

One of us, BRS, would like thank M. Gyulassy for in- 
structive discussions. This work has been supported by 
the U.S. Department of Energy. 
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TABLE I. Properties of the fireballs. 



EOS-I EOS-I EOS-Ib 



Initial parameters 



Rel. fraction of thermal energy 


0.55 


0.60 


0.60 


in the central fireball, Kl 








Longitudinal extension of the 


1.20 


1.20 


1.20 


central fireball, A [fm] 








Rapidity at the edge of the cen- 


1.00 


1.05 


1.00 


tral fireball, va 









Rapidity at maximum of initial 


0.80 


0.85 


0.85 


baryon y distribution, y m 








Width of initial baryon y distri- 


0.32 


0.32 


0.32 


bution, a 








Freeze-out temperature, Tf 


139 


116 


139 


[MeV] 








Output 








Max. initial energy density, ea 


15.3 


16.8 


16.8 


[GeV7/m 3 ] 








Max. initial baryon density, 


3.93 


4.11 


3.90 


nT x [fm- :i ] 








Rel. fraction of baryons in cen- 


0.71 


0.78 


0.73 


tral fireball, f„ B 








Freeze-out energy density, e/ 


0.292 


0.127 


0.127 


\GeVI fm 3 ] 

[l_re V 1 J III J 








Max. transverse velocity at 


0.46 


0.61 


0.58 


freeze-out, v^ ax [c] 








Lifetime of fireball, t ma x [fm/c] 


13.1 


20.2 


20.5 


Lifetime of QGP, tQGP [fm/c] 


2.4 


2.6 


6.0 




£ [GeV/fm 3 ] 

FIG. 1. Ratio of pressure and energy density, P/e, and 
temperature, T, as functions of e, for the equations of state 
EOS-I (solid lines), EOS-II (dashed lines), EOS-III (dotted 
lines), and EOS lb (dashed-dotted lines), respectively. The 
dots in plot (b) correspond for each EOS to the starting values 
of P/e with respect to the achieved initial maximum energy 
density ea at transverse position r± = 0, (dots correspond 
to Tf — 139 MeV whereas the open circle corresponds to 
Tf — 116 MeV). The squares indicate the final values of P/e 
at breakup energy densities, ef. The dots A, B, C in plot (c) 
indicate the relationship between the temperature and the 
energy density at the late hadronic stage of the fireball ex- 
pansion. 
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FIG. 2. (a) Rapidity spectra and (b) transverse mass spec- 
tra, l/m±dN/dm±dy, of negative hadrons, hT , (c) rapid- 
ity spectra of protons (without contributions from A decay) 
and (d) transverse mass spectra, l/m±dN/dm±, of protons 
(including contributions from A decay), p, positive kaons, 
K + , and positive pions, ir + , respectively. The solid (double 
dotted-dashed [dashed-dotted]) lines indicate the results of 
the calculations when using equation of state EOS-I (EOS-I 
[EOS-lb]) with Tf = 139 (116 [139]) MeV. The open circles 
represent preliminary data taken by the NA49 Collaboration, 
whereas the filled squares represent final data taken by the 
NA44 Collaboration. 
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FIG. 3. Isothermes for the relativistic fluids governed by 
EOS-I, and EOS-lb at r± = 0, respectively. For EOS-I the 
lines (beginning with the most outer curves) correspond to 
temperatures, T = 116 MeV, 120 MeV, 140 MeV, 160 MeV 
... etc. For EOS-lb, the outer lines correspond to a tem- 
perature, T = 140 MeV , and each successively smaller curve 
represents a reduction in temperature by AT = 20 MeV . The 
lines c — 1 represent the light cone. 
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FIG. 4. Projections of BEC functions for tt + tt + pairs 
emerging from the low p± (horizontal and vertical) accep- 
tance setting of the NA44 detector [8]. The data points are 
data taken by the NA44 Collaboration [8]. The solid [dashed 
[dotted [dashed-dotted]]] lines correspond to the calculations 
using EOS-I [EOS-II [EOS-III [EOS-lb]]] with T f = 139 MeV, 
and the double dotted-dashed lines correspond to the calcu- 
lation using EOS-I with Tf = 116 MeV. The values of the 
inverse width parameters, Ri (i = l,o,s), have been obtained 
from the Gaussian Bertsch-Pratt parametrization (see text). 
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